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Abstract 

The elastomer properties (large deformations, damping) 
male their use more and more common in industries such as 
aerospace, automotive, construction and civil engineering or 
even the entertainment industry... Modeling the behavior of 
such materials is highly nonlinear, the nonlinearities are 
both geometric (due to large deformations imposed) and 
behavioral (behavior laws used are nonlinear). This paper 
presents a detailed description of the numerical implemen- 
tation of incompressible isotropic hyperelastic behavior. In 
this study, the analysis of large deformation problems of 
Ogden's hyperelastic is based on the finite element method 
(FEM) and mathematical programming. To solve the 
problem of balance, we suggest using the combined 
Newton-Raphson/ Arc-length procedure. A typical example 
is presented to illustrate the performance of this formulation. 
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Introduction 

The theoretical studies of hyperelastic materials have 
mobilized many researchers for decades; see eg the 
work of Treloar [l],Gent [2] and Frakley [3]. However, 
the first important works were proposed by Mooney 
[4] and Rivlin [5]. Valani and al. [6] proposed to write 
the strain energy in a form separated following the 
main directions. This led to the widely used model of 
Ogden [7, 8]. Recently, different models have been 
developed for hyperelastic materials especially in 
biomechanics (Zulliger and al. [9]). 

The objective of this paper is to present the model 
most widely used hyperelastic behavior and in 
particular the isotropic model. To do this, the different 
results of basic mechanics of the major changes are 
outlined in the first place. Are recalled and the various 
tensor quantities used in the modeling of continuous 



medium and medium formulation based on isotropic 
hyperelastic tensor structure. In a second step, the 
most widely used hyperelastic potential is exposed. 

Description of Large Strains 

Elastomers, with hyperelastic behavior, usually work 

in large strains. Note x and X positions vectors of a 
particle P of a deformable body in the current 
configuration (where the solid occupies the volume Q ) 
and initial (where the solid occupies the volume Q ). 

We will use the same reference {o,e, l ,e 2 ,e, i ) for the 
initial configuration Co given at time to and deformed 
configuration G at time t (Fig.l). 




FIG. 1 INITIAL AND DEFORMED CONFIGURATIONS 

For a Lagrangian description, the movement of the 
body can be defined relative to a reference 
configuration Co, by a vector function: 

O: _° ' _ (1) 
[X H> ! = <5>{X,t). 

By introducing the displacement vector u and we 
write the equation (1) in the equivalent form: 

x = x+n(x,t). (2) 

Order to describe the geometric transformations 
associated with these large strains, we introduce the 
deformation gradient tensor F defined by: 
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(3) 



(4) 



d$>.(X,t) 8x. du. 
F..(X,t) = — = — l - = S..+ — -, 

■ jK ' ex. ex. " ax. 

j j j 

Where 

Sy denotes the Kronecker symbol. 
Or equation (3) can be written in matrix form: 
F = I + Vu, 
Where 

I is the unit tensor, 

V is the gradient tensor of displacements. 

Due to the large displacements and large rotations, the 
strain tensor of Green-Lagrange E was adopted to 
describe the nonlinear relationship between the 
deformation and displacement: 



1 



E= l -(C-I), 



(5) 



Where 

C — F T F is the strain tensor of Cauchy-Green right. 

In the case of a hyperelastic law, there is a strain 
energy density W which depends on one of the strain 
tensor and whose derivative gives the second stress 
tensor of Piola-Kirchhoff S: 



dW „ dW 



BE, 



(6) 



In the particular case of isotropic hyperelasticity [10], 
equation (6) can be written as follows: 



S = 2 



Where 





' + 


( dW T dW^ 

+ / i 

v 57, dl 2j 


dw c 


3 dl 3 




dl 2 



(7) 



/ (i = 1, 2, 3) are the invariants of the Cauchy-Green 
tensor on the right, C, such as: 

/, = tr(C) ; I 2 = ^[if -fr(C 2 )] and 7, = det(C), (8) 
Where 

tr(*) and det(«) indicate, respectively, the trace and 
the determinant of a tensor. 

Note that S has no physical meaning but it is 
symmetric and purely Lagrangian, and from which 
one can determine the expression of the Cauchy 
tensor <r : 



det(F) 



-FSF . 



(9) 



For the Ogden model [7] for incompressible materials, 
the energy density is expressed in terms of 
eigenvalues % (/=1, 2, 3) of the Cauchy-Green tensor 

right by: 



w{\,x 2 ,i ) ) = Y 4 —{K k +K l -3), 



(10) 



where the constants are the material parameters and m 
is the total strain energy terms. 

We recall that the incompressibility of the material 
results in det(F) — 1 ou s^A^ — 1. Hence, the energy 

density can be expressed only by \ and A 2 : 



f 



1 



-3 



(11) 



Variational Formulation 

Solving the problem balance of incompressible hyper- 
elastic is reduced to problem of minimizing the 
functional *P : 



Min 



*¥(u,S) = \ n (S.E(u)-f.u)dn-\ r t.udT 
Subject to 



qv "I j r, 

Conditions of incompressibility, 
Boundary conditions (u = u on T u ). 



(12) 



Where 



Q is an isotropic structure, subject to surface 
forces 1 on the surface r t , 

/ is the gravity forces, 

U is displacements imposed on the surface T u . 

The nonlinear equations (12) are solved numerically 
with the finite element method whose the approxi- 
mation of the displacement field is defined by the 
relation: 

u(x) = N(x)U with s = V(N(x))U, (13) 
Where 

x =< x, y, z > are the nodal coordinates, 

N(x) is the shape functions matrix, 

U is the vector of unknown nodal displacements. 

Therefore, the functional *P takes the discretized form 
as follows: 
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W(U) = ^[S.E(N.U)-N T .f.u)dn-^ N T .t.UdT. 



(14) 



SSI 



Thus, the discrete system (14) can be solved iteratively 
by using the combined Newton-Raphson/Arc-length 
procedure (see [11] for details on the implementation 
of this method in the context finite element). 

Numerical Example 

In this example, we consider a plate with a central 
circular hole in the stress plane and tension is sought 
under a constant displacement up to 100% elongation. 
Geometry (L * D = 25.4 x 14.24 cm 2 ) with a unit 
thickness is shown in Fig. 2. For more details see 
reference [12] . Due to symmetry of geometry, only one 
quarter of the plate is modeled, as shown below. 





FIG. 2 GEOMETRY AND LOADING (LEFT) AND BOUNDARY 
CONDITIONS WITH MESHES (528 ELEMENTS T3) (RIGHT) 

UX 
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FIG. 3 DEFORMATION ALONG THE X AXIS (TOP) AND ALONG 
THE YAXB (BELOW). 
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FIG. 4 BOVALUES OF STRESS ALONG THE X AXB (TOP) AND 
ALONG THE Y AXB (BELOW). 

Conclusions 

A computational method by using the mathematical 
programming and the analysis approach hyperelastic 
is developed here to investigate the large deformation. 
A finite element model in conjunction with the 
combined Newton - Raphson / Arc - length procedure 
was adopted allowing investigating interesting 
applications. The obtained results illustrate the 
effectiveness of our approach developed. 

This work can be extended in the future by taking into 
account other parameters as frictional contact with 
heat transfer. 
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